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ABSTRACT 

In the last few years new classes of extremely bright supernovae have been discovered, but 
their rates are so small that models either fail to produce any or dramatically over-produce the 
event rates. These super-luminal supernovae tend to occur almost exclusively is relatively low- 
mass galaxies that undergo active star formation. In the same type of galaxies another high 
energy phenomenon occurs, which are the long-duration gamma ray bursts, which are associated 
with another type of very energetic supernovae, the SN Ic-peculiar. We argue that both the super 
luminal supernovae and the long-duration GRBs are exclusive products of dynamical interactions 
and collisions in young dense star clusters, which are abundant in dwarf galaxies with active star 
formation. We present a model that explains how these different types of explosive events can be 
produced and show that this model can explain their observed rates. In our model the different 
types of super luminal supernovae and the long-duration gamma-ray bursts are related, in being 
a natural consequence of the dynamical evolution of dense star clusters. 

Subject headings: gamma-ray burst: general - supernovae: general - globular clusters: general - galaxies: 
star clusters - galaxies: starburst 



1. Introduction 

In the past decades several new and rare types 
of extremely bright and peculiar supernovae have 
been discovered: (i) The broad-line Type Ic- 
peculiar supernovae that are associated with long- 
duration GRBs ("long GRBs" hereafter). These 
were the first new and ex tremely brigh t type 



of originally very massive stars, and the prototype 
SN1998bw/GRB98042 5 ejected of order half a so 



of supernovae discovered (|Galama et al.l Il998t 
Gehrels fc Meszarod 2012 ) . They have no H and 
He in their spectra and are characterized by ex- 
tremely large outflow velocities (~ 40,000km/s), 
implying very large kinetic energies (~ 10 52 ergs). 
They are related to the explosions of rapidly 
rotating pure C arbon-Oxygen (C O) stars with 
masses > 5M ICano et all l |201lh - bare cores 



lar m ass of Ni 56 (jlwamoto et al.Hl998t ICano et al 
1201 lh . (ii) The so-called Super-luminous Super- 
novae (SLSN), a new class of supernovae dis- 
covered with the recent large-scale surveys for 
transients. The several tens of extremely ener- 
getic and bright SLSNe have bolometric luminosi- 
ties up to some 50 tim es those of type la super- 
novae (|Gal-Yaml 120121 ). There are at least three 
classes of SLSN: the SLSN-I which lack hydrogen 
in their spectra, the SLSN-II which do have H 
in their spectra and the SLSN-R which have a 
long light-curve tail powered by the radioactive 
decay of a large amount (~ 5M Q ) of Ni 56 (see 



1 



Gal-Yaml (|2012h for a review). Both the SLSN 
and the long GRBs with their Type Ic-peculiar 
SNe have in common that: (1) They are extremely 
rare: The rate of long GRBs (lCT 7 /Mpc 3 /yr) is 
some 10 3 times lower than the core-collapse SN 
rate. The combined rate of the SLSN is of or- 
der 10 _8 /Mpc 3 /yr, some 10 4 times lower than the 
core-collapse SN rate, which imply that an ex- 
tremely rare type of stellar evolution is required 
to produce these types of events (see Tab.[l}. 
(2) Both the long GRBs and SLSN occur al- 
most exclusively in small star-forming galaxies 
(SMC/LMC-like or smaller). Only one in a rep- 
resentative sample of 42 long-GRBs is hosted by 
a grand-design spiral galaxy, and all the 41 others 
fall on o ptically bright spots o f small star-forming 
galaxies ( Fruchter et al.ll2006h . Studies of nearby 
small star-forming galaxies show that such bright 
spots are clumps of massive O- and Wolf-Rayet 
(WR) stars. For example NGC 3125 has a num- 
ber of such c lumps that show a mixture of O- and 
WR spectra jHadfield fc Crowther|[200^ . Studies 
of these clumps show that such small galaxies may 
harbor as many as 10 4 O- and WR stars, which 
are concentrated in a small number (3 to 6) young 
massive star clusters. 

2. The connection between super-luminal 
supernovae and long gamma-ray bursts 

The SLSN share the property of the long GRBs 
to occur almost exclusively in small star-forming 
galaxies. The only two SLSN-II that reside in 
larger Milky- Way-type galaxies were fo und very 
close to the nucleus of their hosts and (|Gal-Yam 
2012J) remarks that this "suggests that to pro- 
duce SLSN perhaps special conditions are required 
that are unique to this environment (e.g. circum- 
nuclear star-forming rings), somehow mimicking 
the conditions in star-forming dwarf galaxies." In- 
deed, the inner few hundred parsecs of the bulges 
of most spiral galaxies are undergoing nuclear star- 
bursts. Seyfert galaxies, and e ven the Milky- Wa y 
Galaxy, are prime examples (jConti et al.l 120081) . 
Within 100 pc from the center of the Milky- Way 
Galaxy many dense young massive star clusters 
are present, of which the Arches and Quintu- 
plet clusters are the most p rominent, both host- 
ing > 10 4 stars (|Figerll2004h . Even more impor- 
tant than mass (as we will argue below) is their 
very high stellar density, which in the center of 



Arches exceeds 10 5 M Q /pc 3 . Such a high stel- 
lar density seems to be also a characteristic of 
the young clusters in small star-forming galax- 
ies like the LMC, where for example the central 
part of the R136 clus ter has such a high density 
(jAndersen et al.l l2009h . The reason for the clus- 



ters in small star-forming galaxies to have such 
high stellar densities may be related to the tur- 
bulent velocity structure which leads to a shorter 
free-fall times cale of the gas than in the disks of 
large galaxies (jElmegreen et al l2012h . We suggest 
that the dense torus of in-spiralling gas accumu- 
lating in the central few hundred pc of the bulges 
of spiral galaxies may also have a turbulent veloc- 
ity structure, due to the high local star formation 
and associated high supernova rate. 

It thus appears that both the long GRBs/Type 
Ic-peculiar SNe and the SLSN solely occur in re- 
gions of galaxies where dense young massive star 
clusters are present. This strongly suggests that 
both types of objects are the products of evolu- 
tionary processes that are unique to young mas- 
sive star clusters, and do not occur anywhere else 
(such a suggestion was casually made to one of us 
in 2006 for long-GRBs by S.R. Kulkarni). 

2.1. Conditions required for producing 
long GRBs 

The general consensus is that the condition re- 
quired for producing a long GRB is: the collapse 
of a rapidly rotating practically bare CO-core of 
a massive star. There is strong observational evi- 
dence that the GRB is produced by a narrowly col- 
limated relativi stic jet with a Lorentz-fact or of or- 
der 10 2 to 10 3 (|Gehrels fc Meszarosl[2012l ). There 
are two models fo r producing su ch jets: (i) The 
collapsar model of ( Woosleylll993 ) , in which a very 
massive rapidly rotating core collapses to a black 
hole, but the core has so much angular momen- 
tum that not all of the core matter can at once 
disappear into the black hole, and part of the core 
matter temporarily forms a disk of nuclear matter 
around the black hole. Viscous and/or magnetic 
dissipation on this disk drives a relativistic jet 
which produces the GRB (jMacFadven fc Wooslev 
1999f ): (ii) An alternative collapse model is one 
in which the rapidly rotating core collapses to an 
unusually strongly magnetized neutron star (mag- 
netar) that is spinning with a period of order of a 
millisecond. The spin-down energy loss by mag- 
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netic dipole emission and the relativistic wind of 
such an extreme pulsar is so gigantic that it will 
spin-down to long periods on a timescale of min- 
utes to hours and produce energetic electromag- 
netically powered relativistic jets along the rota- 
tion axis (IWooslevll2010HKasen fe Bildstenl [2O10t 
iGehrels fc Meszarosll2012f )~ 

Both these models require that in order to pro- 
duce a long GRB, the collapsing massive CO 
core must have very hi gh angular momentum 
(jWooslev fe Blooml 120061 ). Two possible ways 
have been suggested for the core of a massive 
star to keep high angular momentum through- 
out its life, either (a): very low metallicity 



(e.g. lYoon fc Langerl (120051)') or (b) evolution in a 
close bi nary system (e.g. van den Heuvel fc Yoonl 
(|2007l ): iFrver et al.l (|2007l )). In the first type of 
models it is argued that low metallicity gives weak 
stellar winds such that the winds do not carry off 
much angular momentum and the star keeps high 
angular momentum throughout its life. Its rapid 
rotation, in these models, keeps the star com- 
pletely mixed, such that it evolves homogeneously 
and may in the end become a rapidly rotating 
CO star that collapses. However, the require- 
ment of very low metallicity is not fulfilled for 
many of the host galaxies of long GRBs, because 
their metallicities tu rn out to range from 0.01 
all the wav to Z^ (IWolf fc Podsiadlowskil [2007 



IGehrels fc Meszarosl 2012t ). Therefore low metal- 
licity cannot be the defining cause of the occur- 
rence of long GRBs, since whether or not a long 
GRB will occur in a galaxy appears not to be 
strongly dependent on its metallicity. 

In close binary models tidal f orces keep the star 

in syn chronous (rapid) rotation ([van den Heuvel fc Yoon 
2007), or a rapidly rotating m erger product is pr o 



duced (Helium merger GRB; IFrver et al.1 (120071 )1. 
The first model has great difficulty reaching the 
specific angular momentum required for making 
a long GR B, and may only do so i n very special 
cases (e.g. Ivan den Heuvel fc Yoonl (|2007l) ). Since 
massive binaries are found throughout the disks 
of spiral galaxies, both models, even if they could 
succeed, would be expected to produce many long 
GRBs in disks of spirals, contrary to the observa- 
tions. 

While normal close binary evolution, including 
tidal synchronization, may not be able to reach 
and maintain the specific angular momentum re- 



quired for producing a GRB, a CO-core produced 
in an off-center collision of two evolved massive 
stars in a dense cluster, may have very high angu- 
lar momentum, sufficient to produce a GRB. 

3. The link with dense young star clusters 

The majority of massive stars form in dense 
clusters. The most massive star m max in a clus- 
ter of mass M r1 is m max ~ 1.2(M c i/M Q ) ' 45 M m 



(jWeidner fc Kroupal 12004) with a maximum of 
about 150 M Q ( Figerl 2005). Such a star sinks to 
the cluster center by dynamical friction in a small 
fraction ((m)/m max ) of the two-body relaxation 
timescale i r i x , although the star cannot reach the 
cluster center within it's own orbital time scale in 
the potential of the cluster. 

Massive stars (of mzAMS ~ 25 M Q ) develop 
massive CO-cores (m co > 4M ) near the end of 
their fuel processing lifetime; we calculate the mo- 
ment until the CO core develops by performing a 
series of si mulations using the MESA stellar evolu- 
tion code dPaxton et al.ll201ll) within th e AMUSE 
framework ( Portegies Zwart et al. 20131 ). and this 
fits 



tr 



52(m ZA Ms/M ) 



-0.635 



[Myr]. (1) 



For a sufficiently dense cluster, t co of the most 
massive star is longer than the dynamical fric- 
tion time scale, and the star develops a CO core 
after it arrives in the cluster center. Many of 
these massive stars are born in close binary sys- 
tems (|Sana et al.l |2012| ). the majority will evolve 
according to regular Roche-lobe overflow. A sin- 
gle massive star, or one in a binary system, will 
upon arrival in the central portion of the star 
cluster, acquire a companion to for m a binary 
or triple system ( Heggie et al. 19961) . A newly 
formed binary will at first be rather wide, with 
a binding energy comparable to the mean ki- 
netic energy of the stars, or ~ lkT. Repeated 
interactions with other cluster members drive 
the harde ning of the binary to ^ 100 kT in 
< 0.2i rlx (|Portegies Zwart fc McMillan! I2002T) . A 
side effect of this hardening process, mediated 
by exchange interactions, is the ejection of mas- 
sive stars with high velocity as OB runaways 
(|Fujii fc Portegies Zwartl 120111 ) . This hardening 
process eventually causes the two most massive 
stars in the cluster to occupy the same binary 
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(jGaburov et al.l 120081 ). although more bodies are 
involved in the process ([Tanikawa et al.ll2012l ). 

The chance of a traffic accident during the quick 
dynamical evolution of this central sub-system is 
very large, in particular if the configuration resem- 
bles a hierarchical system in which the two most 
massive stars are in a relatively tight binary that is 
orbited by a less massive outer star. Such systems 
with inclination i £ 50° are kn own to be subject 



to the Kozai effect (]Kozail ll962), in which the ex- 
change of angular momentum between an inclined 
outer orbit and the tight inner orbit drives the 
latter to extremely high ( £ 0. 9) eccentricity. The 
time scale for the Kozai cycles (jKinoshita fc Nakai 
1999) is much shorter than the core relaxation 
time of the cluster, and the two inner stars even- 
tually experience a collision. 

The consequence of a collision depends on the 
evolutionary states of both stars. If sufficiently 
close in mass, the two stars could both have si- 
multaneously well developed CO cores. A dynam- 
ical (or Kozai-driven) off-center collision between 
these stars results in common-envelope evolution 
in which the two CO-cores spiral towards each 
other and the collision energy combined with the 
released orbital binding energy drive off the El- 
and He- envelope. Mergers of the two inspiralling 
CO cores produces a critically rotating massive 
CO-st ar which resembles the progenitor of a long 
GRB ( Cano et al.l |20 1 lh . During the explosion 
of the £ 5 M Q rapidly rotating helium depleted 
merger product, a considerable amount of Ni 56 is 
expected to be ejected in a Type Ic-peculiar su- 
pernova (jCano et al.ll201ll ). 

The time interval in which the star has a CO 
core, however, is extremely short, and at first 
glance it seems unlikely to have both stars of an 
interacting binary in this state. The cluster, how- 
ever, assists here in bringing the two most mas- 
sive stars together in a binary star. The life- 
time of the CO core (m co £ 4M ) of a massive 
m zams > 25 M star is 4000 to 14000 years. We 
constrain the duration between the development 
of the CO-core and the moment of the supernova 
with the earlier mentioned simulations, and find it 
to fit 

dt co ~ 0.15(m zams /M Q )-°- 70 [Myr]. (2) 

Inverting Eq.[2] gives the allotted range in 
masses between primary and secondary star for 



which both have a well developed CO core at the 
same tim e. By integrating the initi al stellar mass 
function (IKroupa k, Weidnerl 120031 ) we calculate 
the probability that two stars in a cluster of mass 
M c \ both have a massive CO core. The integration 
was performed b y adopting a Schechter function 
(|Schechterlll976l ) for clusters masses with a mini- 
mum mass of M c \ = 10 3 M© and a characteristic 
mass of I0 6 M©. To bracket the event rate for the 
Milky- Way Galaxy we also integrate over a distri- 
bution of cluster radii. For the size distribution 
of the clusters we fitted the observed distribution 
of cluster sizes (taken from Tabs 2, 3 and 4 of 
Portegies Zwart et al.l (|2010h ) to a log-normal dis- 
tribution, which gave a satisfactory fit for a mean 
radius of 5 pc and a dispersion of 3 pc. 

4. Event rates for SLSN and GRBs in 
young star clusters 

Normalized to the core-collapse supernova rate 
(derived from counting the number of stars be- 
tween 8 and 25 M©) we then obtain a relative fre- 
quency of 1Zi c - p ~ 1.1 • 10~ 4 , which is very low 
compared to the obs erved rate of 7£i c _ D ~ 10" 



(jGalama et al.ll2000h . The rate in blue-compact 



dwarf galaxies may be higher due to the over abun- 
dance of relatively low-density star clusters in the 
Milky- Way Galaxy, which we used to constrain the 
size distribution of the star clusters in our simu- 
lations. Reducing cluster radii to a mean of 2 pc 
but keeping the same dispersion in the log-normal 
distribution results in a rate of 7?.i c _ p (0.14 — 
14) • 10~ 3 , which brackets the observed rate. Here 
the strict upper limit is calculated by allowing also 
less massive stellar pairs (but with mco > 4M©) 
to contribute to the supernova type Ic-peculiar 
rate. This derived rate brackets the observed rate 
and is consistent with the observed rate for long 
GRB (jGal-Yamllioil . We summarize the ob- 
served and our derived rates in Table[TJ With 
the same integration we calculated that stars with 
wzams ~ 25 M© collapse in a type Ibc supernova 
at a rate of IZihc — 0.25 ± 0.01, which is consis- 
tent with the observed rate of IZibc — 0.22 ± 0.06 
(jCappellaro et al.lll999l ). 

The above scenario does not work for clus- 
ters that experience core collapse before the 
first supernova, because this will initiate a col- 
lision runaway producing a supermassive star 
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(jPortegies Zwart k McMillaij|2002l ). This process 
is arrested by stellar-wind mass loss and the first 
supernova. The amount of mass that is accumu- 
lated in the collision runaway product depends on 
the time of core collapse co mpared to the lifetime 



of the most massive star (jPortegies Zwart et al 



20061 ). As a consequence, the event rate depends 
on cluster mass and size. 

The collision rate during the time between core 
collapse and the supernova explosion of the colli- 
sion runaway produ ct determines the maximum 
mass of the latter (jPortegies Zwart fc McMillan! 
2002). In relatively low-mass star clusters M c \ ~ 
10,000M o to 30,000M Q the collision runaway 
product can reach a mass of 150 M^ to 260 M Q 
(jPortegies Zwart fc van den Heuvelll2007j ). These 
stars coll apse in a luminous electron-capture su - 
pernova (jScannapiecol 120091 : ICooke et all l2012h . 
giving rise to a SLSN-R, as these supernovae pro- 
duce large a mounts of Ni 56 , li ke was proposed for 
SN2007bi bv lPan et all (j2012f ). 

Only the densest (virial radius ^ 0.3 pc) star 
clusters contribute to this rate. An environment 
with a larger proportion of dense and massive 
clusters, like in blue-compact dwarf galaxies, will 
have a higher proportion of SLSN-R compared to 
core collapse supernovae. For the cluster size dis- 
tribution as observed in the Milky- Way Galaxy 
we calculated a rate of 7^slsn-r ~ 2.3 • 10~ 6 , 
which is much lower than the observed rate of 
^slsn-r = 2 • 10~ 5 . The rate we derive for 
the population of blue-compact dwarf galaxies is 
7^slsn-r — 2.1 • 10 -5 , which is close to the ob- 
served rate. 

In clusters of > 30,000M Q , if sufficiently 
dense the collision runaway can grow to a mass 
> 260 Mq. We suggest that these "super stars" 
produce SLSN-I/II by collapsing to a bla ck hole 
of intermediate mass IScannapieco ( 2009() . The 
mass of the collision runaway can reach values of 
up to a few 10 3 M (jPortegies Zwart fc McMillan 
2002). By the time the star experiences a su- 
pernova it has she d most of its mass aga in in a 
dense stellar wind ( Yungelson et al.ll2008l) and it 
is uncertain how much mass eventually collapses 
to the black hole. Integrating over the mass and 
size distributions for star clusters we derive a rate 
of ftsLSN-i/n ^ 6.1 • 10~ 5 for the Milky-Way 
population. By adopting the same size distri- 
bution and mass distribution of star clusters as 



we did before for the population of clusters in 
blue-compact dwarf galaxies we arrive at a rate of 
^•slsn-i/ii — 1-8 • 10~ 4 , which again is similar to 
the observed rate for combined types SLSN-I and 
SLSN-II. 

We do not make the distinction between type 
SLSN I and type II here but derive the total 
rate. Upon each collision several M Q of hydro- 
gen is injected into the collision runaway, but 
this mass is blown away in the copious stellar 
wind in a few 10 4 years. A collision between 
the runaway and a hydrogen rich star shortly be- 
fore the supernova of the former was prop osed by 
( Portegies Zwart fc van den Heuvel 120071 ) to ex- 
plain the SLSN-II 2006gy wh ich occurred cl ose to 
the nucleus of a large galaxy ( Quimbv 20061) . The 
ratio between the timescale on which fresh hydro- 
gen is injected into the collision runaway and the 
time required to deplete the newly acquired hydro- 
gen envelope determines the ratio of SLSN type II 
relative to SLSN type Is. The observed compa- 
rable rates of type SLSN-II relative to SLSN-I is 
consistent with this regime of colli sional growth 
(jPortegies Zwart fc McMillan! I2002I ). Both the 
rates for SLSN-R and for SLSN-I/II increase if a 
higher proportion of star clusters are born with 
high density, like is the case for clusters in blue- 
compact dwarf galaxies compared to the Milky- 
Way. 

In the collapsed core of a dense young star 
cluster an intermediate mass a black hole is likely 
to be accompanied by another star, or otherwise 
it will acquire one within a core relaxation time 
scale. The orbital period of such a b inary typically 



is in the range of 50 to 500 days (jPatruno et al 
20061) . Once the captured companion over- 
fills its Roche-lobe and transfers mass to the 
intermediate-mass black hole the binary becomes 
visible as an ultra luminous x-ray source, much 



like t he observed syst ems M82 X-l (IKaaret et al 
200 ll ). NGC1 313 X-2 dZampieri fc Patrunoll201ll) 



and HLX-1 dWebb et al l |2012|), NG C5408 X-l 



(IStrohmaverll2009l) and NGC7479 X-l (jVoss et al 
l201lj ). The observed periodicity in M82 X-l 
(62 days), NGC5408 X-l (115days) and HLX-1 
(388days) and their x-ray fluxes are consistent 
with a cluster member being captured by an in- 
termediate mass black hole and feeding the latter 
via a dense stellar w ind or Roche-lobe overflow 
( Patruno et al. 20061 ). We expect such objects to 
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reside in clusters of > 30, 000 M Q . 
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Table 1: Event rates for the families of supernovae. Observed rates (second column) are from (|Gal-Yamll2012l) . 
The rates are normalized to the core collapse supernova rate, which for the simulations was calculated sclf- 
consistently. The best values are from our adopted Schechter mass function with a characteristic mass of 
10 6 M© and with a log-normal size distribution with mean of 2 pc, which represents the star clusters in 
blue-compact dwarf galaxies. The rates for the Milky- Way are calculated for a sample of young star clusters 
consistent with those observed in the local group. 







model 


model 


SN type 


observed 


best value 


Milky-Way 


SN Ibc 


0.22 


0.25 


0.25 


SN Ic-peculiar 


1 ■ 10" 3 


(0.14-14) ■ 10~ 3 


1.1 • 10~ 4 


SLSN-I/II 


1.7- 10~ 4 


1.8 • 10~ 4 


6.1 • 10~ 5 


SLSN-R 


2 • 10" 5 


2.1 • 10~ 5 


2.3- 10~ 6 
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